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We demonstrate that whispering gallery mode resonators can be utilized to generate optical 
frequency combs based on four wave mixing process at virtually any frequency that lies in the 
transparency window of the resonator host material. We show theoretically how the morphology of 
the resonator can be engineered to produce a family of spectrally equidistant modes with anomalous 
group velocity dispersion appropriate for the comb generation. We present experimental results for 
a frequency comb centered at 794 nm to support our theoretical findings. 



Since their introduction over a decade ago, optical 
frequency combs have revolutionized many applications 
across the fields of science and metrology, and have 
opened new vistas in sensing and spectroscopy. The 
conventional approach for optical comb generation with 
femtosecond lasers was recently augmented with optical 
combs generated using whispering gallery mode (WGM) 
resonators made of materials possessing Kerr nonlinear- 
ity [1-12]. These small structures address one of the 
shortfalls of femtosecond laser frequency combs, namely 
their size and relatively complex architectures that limits 
their range of applicability outside research laboratories, 
and thus restricts their applications. The second shortfall 
of optical frequency combs related to the frequency of op- 
eration and has not been effectively addressed. Though a 
number of optical combs covering selected spectral range 
from UV to IR have been demonstrated, there is no uni- 
versal approach that can provide an optical frequency 
comb centered at any desired wavelength. We propose a 
solution of the problem in this paper. 

The generation of a phase locked frequency comb with 
an optically pumped WGM microresonator is mediated 
by hyper-parametric oscillation (modulational instability 
[16-18]) resulting from four wave mixing (FWM), and 
by group velocity dispersion. In the FWM process, two 
pump photons and two sideband photons are involved. 
The generation of two sidebands symmetrically around 
the excited mode starts when the power of a continuous 
wave laser pumping the mode exceeds some threshold 
value [15]. Increasing the pump power results in gen- 
eration of multiple sidebands via the hyper-parametric 
process, producing a phase locked group of tines in the 
optical comb. Depending on the dispersion of the WGM 
spectrum, the initial generation of sidebands separated 
by multiple free spectral range (FSR) could have a lower 
threshold compared with that of closely separated side- 
bands [3, 12]. 

Group velocity dispersion (GVD) also plays a critical 
role in the generation of frequency combs. To explain 
this, let us recall that the threshold condition for the 
average number of photons, N, in the pumping mode 
necessary to start the FWM oscillation process, which 
generates the Kerr frequency comb in a WGM resonator, 



reads as [15] 

4(g 2 N 2 - 7 2 ) = [2(wo - u) + (w+ + U- - 2oj ) - 4<?A>] 2 , 

(1) 

where 7 is the full width at the half maximum of the 
optical modes (we assume that the modes are identical 
in loss), wo and ui± are the frequencies of the optically 
pumped mode and the first pair of the modes where 
FWM sidebands are generated at, to is the carrier fre- 
quency of the external monochromatic pump, g is the 
nonlinearity parameter 
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c is the speed of light in the vacuum, V is the mode 
volume, no and H2 are the linear and nonlinear refractive 
indices of the resonator's host material. 

In accordance with Eq. (1) oscillation is allowed if 
at least one of the following two conditions is satisfied: 
i) the optical pump is detuned to the lower frequency 
than the corresponding optical resonance (red detuned), 
luq > ui; and ii) the resonator modes are characterized 
with anomalous dispersion, f} 2 < (2w — w + — w_ ~ 
cf32Ujp SR /no, where lufsr is the free spectral range of 
the resonator). The condition (ii) is usually stronger 
than (i). Even though oscillation can occur in a res- 
onator possessing arbitrary frequency dispersion in the 
mode family participating in the comb generation [16— 
18], it generally starts when the dispersion is small 
and anomalous. Resonators with normal dispersion ex- 
hibit much weaker FWM process, which frequently com- 
petes with stimulated Raman scattering (SRS), which 
has a lower threshold for oscillation [19]. The red de- 
tuning of the pump light with respect to the corre- 
sponding resonator mode is required for phase match- 
ing and compensation of the dispersion, though it in- 
creases the oscillation threshold power, so the SRS pro- 
cess starts first. By contrast, anomalous dispersion sim- 
plifies the phase matching of the FWM process [5] and 
removes the burden away form laser detuning. These 
assertions have been confirmed experimentally. For in- 
stance, for combs observed with WGM resonators, a 
value of (2wo — w+ — w_)/(27r) — —2.6 MHz was mea- 
sured for a 75 /xm fused silica microresonator at 1,550 nm 
[1]. The dispersion of a 0.7 cm CaF2 resonator was 
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(2w -U3 + - w_)/(2tt) ~ -300 Hz [20], and a 0.255 cm 
CaF 2 resonator had (2lu — u + — cj_)/(27t) ~ —680 Hz 
[3]. 

The dispersion fa includes both geometrical and mate- 
rial parts. The geometrical dispersion is usually normal. 
To obtain mode spectra characterized with small anoma- 
lous dispersion one needs to produce a resonator out of 
a material with anomalous dispersion for the given laser 
wavelength. The shape and size of the resonator must 
be properly selected to have an overall small and anoma- 
lous dispersion (\2lu — u> + — uj_\ <C 7). These require- 
ments prohibit generation of Kerr frequency combs with 
an arbitrary repetition rate at an arbitrary wavelength, 
since the repetition rate is determined by the resonator 
size, and the wavelength is determined by the material 
dispersion. For instance, a fused silica resonator is not 
expected to generate Kerr frequency comb if pumped at 
a wavelenght shorter than 1.3 fim. 

Now, consider a spheroidal resonator with two nearly 
equidistant mode sequences characterized with frequency 
intervals (free spectral range, FSR) FSRi — c/(2Trna) 
and FSRp = FSRi(a — b)/b, where I and p are the az- 
imuthal and the vertical numbers of the mode, c is the 
speed of light in the vacuum, n is the refractive index 
of the material, and a and b are the semi-axes of the 
spheroid [14]. Modes belonging to a basic sequence char- 
acterized with FSRi are generally utilized for generation 
of Kerr combs. However, for the purpose of comb gener- 
ation at any arbitrary frequency within the transparency 
window of the resonator material, we propose to use the 
vertical mode family characterized with FSR P , which de- 
pending on the ratio a/b can have an FSR significantly 
different from FSRi ■ An important feature of a vertical 
mode family is that the mode spacing strongly depends 
on the shape of the spheroid, and the corresponding GVD 
can be anamoulous. 



frequency spectrum of the cylinder is given by 
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FIG. 1: Illustration of the cross sections of a spheroid and a 
cylinder, and a WGM resonator being a convolution of those 
shapes. A spheroidal WGM resonator has normal geomet- 
rical group velocity dispersion in vertical modes, and nearly 
equidistant spectra. A cylindrical resonator has anomalous 
group velocity dispersion for its modes but not equidistant 
spectra. The resonator designed as a convolution of those 
two shapes has a nearly equidistant spectrum and anoma- 
lous group velocity dispersion necessary for achieving phase 
matching in the FWM process for any resonator host mate- 
rial. 



The morphological dependence of the dispersion can 
be understood as follows. Consider the convolution of 
a cylindrical and a spheroidal resonator (Fig. 1). The 
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where ki^ q a ~ I + a q (l/2) 1 ' 3 is the azimuthal wave num- 
ber given by the spherical Bessel functions, q is the radial 
wave number of the mode, a q is q th root of Airy func- 
tion (Ai(—a q ) = 0), and L is the hight of the cylinder. 
In the case of I 3> p the spectrum is significantly non- 
equidistant with p, but the group velocity dispersion is 
anomalous 2ui p j tq — ui p +u tq — ui p -i t i q < 0. The nearly 
equidistant frequency spectrum of a spheroid can be ap- 
proximated by [14] 
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so that the dispersion is normal, 2ui Pi i iq — Wp+i j lQ — 
^p— i,Z,a > 0. A resonator being a convolution of the 
cylinder and spheroid has a frequency spectrum given by 
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(a-fc)(l + 2p) 2 
kLq+ 2ab +CP 
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where £ is a numeric parameter depending on L, a, and 
b. The parameter £ can be positive in the case of compa- 
rably small L < 2b, which means that the dispersion of 
the vertical resonator modes is anomalous. In this way, 
the anomalous dispersion can be used to compensate for 
the normal material dispersion in the Kerr comb. 

There are no measurements of the group velocity dis- 
persion for the majority of optical materials, and rigorous 
calculations confirming these prediction are rather math- 
ematically involved. So we performed an experiment to 
verify these estimates. We designed and fabricated a 
truncated spheroidal crystalline CaF2 WGM resonator 
and pumped it with 794 nm light emitted by a semi- 
conductor laser that was self-injection locked [13] to a 
selected WGM. The light was sent in and retrieved out 
of the resonator using a glass coupling prism. The exit- 
ing light was collimated and injected into a single mode 
optical fiber to be further analyzed. The spheroidal res- 
onator with 3 mm and 3/4 mm semi-axes (with a/b « 4) 
was fabricated out of a z-cut CaF 2 wafer by mechanical 
polishing. The loaded quality (Q-) factor of the resonator 
exceeded 3 x 10 9 (126 kHz full width at the half maximum 
of the corresponding WGMs). 

We observed the onset of hyperparametric oscillation 
with 0.1 mW of optical power. Increasing the power to 
2 mW resulted in the formation of a well pronounced op- 
tical frequency comb shown in Fig. (2). The frequency 
envelope of the comb suggests that the vertical mode 
family is involved in this process. The asymmetry of 
the spectrum results from the significant spatial spread 
of the light emitted from the resonator that is selec- 
tively collected with the fiber. Since the group velocity 
dispersion of the material is normal at this wavelength 
(fa = 28 ps 2 /km) and the difference between two consec- 
utive FSR's is positive and large (2uq — u + — uj-)/(2ir) ~ 



3 



21 kHz, but the comb is still generated, we conclude that 
the anomalous dispersion of the resonator compensated 
the material dispersion. The experiment proves our the- 
oretical prediction. 
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FIG. 2: The output spectra of a critically coupled CaF2 res- 
onator pumped with 2 mW of 794 nm light. No stimulated 
Raman scattering is observed. The resonator has 8 GHz 
free spectral range, while the optical frequency comb has 
23.78 GHz repetition rate. The comb harmonics are phase 
locked, so they generate a coherent RF signal (inset) on a 
fast photodiode. 

Note that the repetition frequency of the observed 
comb is approximately three times as large as the ba- 
sic mode frequency FSR. In our case a/b w 4, so that 
FSRp m SFSRi, which is close to experimentally mea- 
sured value FSR p /FSRi = 23.78 GHz/8 GHz = 2.96, 
where 23.78 GHz is the repetition frequency of the comb 
we observed. This is further proof that the frequency 
components of the Kerr comb reported here are not gen- 
erated in the same family of fundamental modes with dif- 
ferent orbital momenta. Rather, the components belong 
to the family of modes with different vertical quantum 
number p and the same orbital momentum I. 

The maximum number of Stokes and anti-Stokes com- 
ponents of a conventional Kerr comb generated in the 
basic mode family [1-12] is limited by the dispersion and 
nonlinearity of the resonator and can be very large. On 
the contrary, the number of the Stokes comb components 
observed in our system is limited by the vertical index p 
of the mode that is selected for optical excitation. In- 
deed, figure (2) represents an optical comb produced by 
pumping the mode with p = 42. The saturated comb 
has a flat top and a sharp edge at the lower frequency 
side, because there are simply no modes with negative 
transverse index. This type of Kerr comb affords much 
more flexibility compared with a comb generated at the 



basic mode sequence: its repetition frequency is large in 
an oblate spheroid, as observed in our experiment, and 
a nearly spherical resonator would generate a comb with 
repetition frequency much lower than the azimuthal FSR 
of the resonator. In general, the repetition frequency of 
a transverse comb depends symmetrically on both the 
resonator semi-axes a and b 
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This property results in the ability for generation of low 
frequency combs even in very small microresonators of 
very low mode volume, and therefore, very low threshold. 
It also ultimately leads to a new kind of RF photonic 
oscillator [21] with extremely low phase noise. 

The asymmetry of the observed comb envelop requires 
a special explanation. A WGM of vertical index p is de- 
scribed by an eigenfunction which has p + 1 maxima in 
the electric field along the vertical axis of the resonator 
[14] . The overlap of the evanescent field of the mode on 
the coupling prism determines the emission pattern. The 
emission pattern of the transverse mode consists of two 
symmetrical lobes with respect to the horizontal plane. 
Since the angular space between the lobes and the fre- 
quency of the comb components grows with index p, com- 
ponents of lower index p are concentrated closer to the 
equatorial plane of the resonator. The optical fiber cou- 
pler that we use is aligned in the equatorial plane of the 
resonator and thus has much higher efficiency for col- 
lecting the modes with smaller p and, therefore, lower 
frequency. This is why the comb as observed has the ap- 
pearance of being asymmetrical. We note that the spatial 
pattern of the comb emission provides a simple way to 
physically separate strongly correlated Stokes and anti- 
Stokes components of the comb for applications where 
correlated photons are desired. 

In conclusion, we have proposed and experimentally 
validated a universal method for the generation of Kerr 
frequency combs at any desirable wavelength, and with 
any desirable repetition rate. The method is based on 
the management of the geometrical group velocity disper- 
sion of whispering gallery modes participating in the pro- 
cess through the proper engineering of the shape of the 
resonator. We have found that the vertical sequence of 
the modes belonging to a specially shaped resonator has 
nearly equidistant spectrum characterized with anoma- 
lous group velocity dispersion. Using these modes we 
were able to generate a frequency comb at 794 nm center 
wavelength in a properly designed calcium fluoride res- 
onator. The conventional Kerr comb produced with the 
basic sequence of the azimuthal modes was absent in our 
observations. 
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